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' We consider tiie quark-meson-model in a finite three-dimensional volume using the 

o ; 

. Schwinger proper-time renormalization group. We derive and solve the flow equa- 

CN ■ 

5h ' tions for finite volume in local potential approximation. In order to break chiral 
! 

^ ' symmetry in the finite volume, we introduce a small current quark mass. The corre- 

. sponding effective meson potential breaks chiral 0(4) symmetry explicitly, depending 

^ \ on a and vr fields separately. We calculate the volume dependence of the pion mass 

\^ ' and of the pion decay constant with the renormalization group flow equations and 

■ compare with recent results from chiral perturbation theory in a finite volume. 
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O ; I. INTRODUCTION 

^ ■ 

^ . The study of QCD in a finite volume has been of interest for quite some time. Accurate re- 

rH ' suits of lattice simulations with dynamical fermions necessitate understanding finite volume 

■ effects. A variety of different methods has been proposed, cf. refs. 
to extrapolate reliably from finite lattice volumes to the infinite volume. Finite volume 

heir own right, because they allow 
[l2 , Q| . The low energy behavior of 



X 
J3 



partition functions for QCD have attracted interest in t 



an exact description of QCD at low energies 
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QCD is determined by spontaneous chiral symmetry breaking [1^, which, however, does not 
occur in a finite volume. If the current quark mass is set equal to zero, in a finite volume the 
expectation value for the order parameter of chiral symmetry breaking vanishes, remaining 
zero even for arbitrary large volumes. The order parameter has a finite expectation value 
only when the infinite volume limit is taken before the quark mass is set to zero. 

The box size L, the pion mass m^r and the pion decay constant /^r are the relevant 
scales for the transition between the regimes with a strongly broken and with an effectively 
restored chiral symmetry As a measure of explicit symmetry breaking, the pion mass is 
of particular importance. It is primarily the dimensionless product m-j^L that determines in 
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which regime the system exists for a given pion mass and volume. In order to study chiral 
symmetry breaking in a finite volume, it is essential to introduce a finite quark mass as a 
parameter that explicitly breaks the chiral symmetry. Such an explicit symmetry breaking 
is quite natural in theories which involve effective chiral Lagrangians. 

QCD at low energy can be studied by a wide variety of approaches which in essence all 
rely on the same fact: Due to spontaneous breaking of chiral symmetry, low-energy QCD 
is dominated by massless Goldstone bosons associated with the broken symmetry. Since 
these Goldstone bosons interact only weakly, the low-energy limit of QCD can be described 
in terms of an effective theory of these fields. A description in terms of effective chiral 
Lagrangians becomes even better if one considers the partition functions in finite Euclidean 
volume. Compared to the light degrees of freedom, contributions of heavier particles are 
suppressed by e~*^^, where M is the typical separation of the hadronic mass scale from the 
Goldstone masses. This separation of mass scales is at the origin of the description of QCD 
with effective theories in terms of the light degrees of freedom only. 

Groundbreaking work has been done by Gasser and Leutwyler Q| in chiral 

perturbation theory, and by Leutwyler and Sr nilg a for the eigenvalue spectrum of the QCD 



Dirac operator 13]. Random matrix theory 
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redicts analytically the volume and 



predict 

HQ 



and the eigenvalue spectrum 
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21[ . Such analytic 



quark mass dependence of the chiral condensate 
which has been well confirmed by numerous lattice results, see eg. 
predictions have been extremely useful as a check for calculations in lattice gauge theory. 

In the context of the renormalization group fRG) approach, most calculations so far 
have been done in a chirally symmetric formulation (221 . l23l . vlA l25| . In such a framework, 
explicit symmetry breaking terms do not affect the renormalization group flow. This means 
that pion contributions to the renormalization group flow come from exactly massless pions. 
Small quark masses can be converted into a small linear term in the meson potential which 
is not renormalized and can therefore be added at the end of the evolution. In this paper, we 
will present a formulation which includes explicit symmetry breaking in the renormalization 
group flow and which can be used also for a flnite volume. 

We use the chiral quark-meson-model, which contains a scalar sigma meson and quarks 
in addition to the pion degrees of freedom. It is well known that the linear sigma model 
alone is not compatible with the low energy vrTr-scattering data ^1^. Due to the presence 
of quarks, the low energy constants of chiral perturbation theory are reproduced [3, ■ 
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The quark-meson-model is evolved with renormalization group flow equations which connect 
different scales. As we will demonstrate, the inverse box size 1/L acts as an effective cutoff 
scale which freezes the evolution when the evolution parameter k < n/L. In the same way 
as the renormalization group flow equations describe the dependence of the results on the 
renormalization cutoff scale k, they also describe the dependence on the additional scale 
imposed by the finite volume. The summation of higher loop graphs in this approach does 
allow to extend the calculation to smaller quark masses and volumes, where perturbative 
calculations lack convergence. Deviations of the calculated pion mass correction in the finite 
volume from chiral perturbation theory are found in this region. Our results agree, however, 
in the case of large pion masses and large volumes. 

The paper is organized as follows: In section m we review recent results from the appli- 
cation of chiral perturbation theory to the problem of finite volume effects. In section IIIH 
we show how the evolution equation for the effective meson potential is modified in a finite 
volume. Details of the numerical evaluation are discussed in section IIVI In section E] we 
present our final results, which are discussed in section IVIl 



II. FINITE VOLUME EFFECTS IN CHIRAL PERTURBATION THEORY 



For finite volume, massless Goldstone bosons dominate the action of a theory with bro- 
ken chiral symmetry. In chiral perturbation theory (chPT) |l4 |, t he pion mass, the pion 
decay constant and the chiral condensate have been calculated |lo[ 11, Q|. The expansion 
parameters are the magnitude of the three-momentum |p| and the mass of the pion as 
the lightest degree of freedom compared with the chiral symmetry breaking scale 47r/^. 

Depending on the size L of the volume and the pion mass m.,r, chiral perturbation theory 
distinguishes between two different power counting schemes. If the size of the box is much 
larger than the Compton wavelength of the pion L ^ l/rriT^, the lowest nonzero pion mo- 
mentum is smaller than the pion mass (pmin ~ ^ <^ ^n) and the normal power counting 
scheme applies ( "p-regime" ) . In this case, the pions are constrained very little by the pres- 
ence of the box, and finite size effects are comparatively small If, on the other hand, 
the size of the box is smaller than the Compton wavelength of the pion, the normal chiral 
expansion breaks down, since the smallest momentum Pmin ~ T" ^ ''^■^ '^'-'^ much larger 
than the pion mass ( "e- regime" ) . In this case, the partition function is dominated by the 
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zero modes. After solving the zero-momentum sector of the theory exactly, one expands the 
finite momentum modes to one-loop order. 

A ver y us eful tool to study the effects of a finite volume on the mass of the pion is Liischer's 
formula 28j. It relates the leading finite volume corrections for the pion mass in Euclidean 
volume to the vrTr-scattering amplitude in infinite volume. Corrections to the leading order 
behavior drop at the least as 0{e~^^) where ffi > ^J'ijlmT^. For the particular case of the 
pion mass, the formula for the relative deviation R\m.^(Ly\ of the pion mass rriT^iV) in the 
finite volume from the pion mass in the infinite volume 7717^(00 ) reads as follows: 

m-j.^L) — m,r(oo) 



i?K(L)] 



m^(oo 

3 



1 1 /"°° / 

dyF{iy)e-V"''^y'^ + 0(6-'^^). (1) 



F is the forward vrvr-scattering amplitude as a function of the energy variable s continued 
to complex values. New results 0, |^ have recently been obtained by combining Liischer's 
formula with a calculation of the scattering amplitude in chiral perturbation theory. A 
next-to-next-to leading order calculation of F alone does not seem to give a reliable and 
satisfactory result. A one-loop calculation using Liischer's formula gives a shift in the pion 
mass, for example, which is substantially lower than the one expected form the full one- 
loop calculation in chiral perturbation theory as performed by Gasser and Leutwyler 
This estimate of the finite volume effects can be improved, if one uses the mass correction 
obtained from Liischer's formula with a Trvr-scattering amplitude including higher orders to 
correct the full one-loop chiral perturbation theory result ^. Since Liischer's formula has 



sub-leading corrections O yexp[— ^y3J2'm^^L]j , the corrections to the leading result increase 
for decreasing pion mass at fixed volume size L with rriTrL. As pointed out by the authors 
in [^, the Liischer formula becomes a less reliable approximation exactly for those values of 
the pion mass for which the chiral expansion converges especially well. 

The renormalization group flow equations do not rely on either box size or pion mass as 
an expansion parameter, and do not require to distinguish between two different regimes. 
They remain valid as long as the lowest momenta and the masses of the heaviest particles 
remain below the ultraviolet cutoff scale Ajjv ~ L5 GeV. The beauty of the renormalization 
group method is precisely that the flow equations connect different scales. In the same way 
as the renormalization group flow equations describe the dependence of the results on the 
infrared cutoff scale k, they also describe the dependence on the additional scale imposed 
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by the finite volume. 



III. RENORMALIZATION GROUP FLOW EQUATIONS FOR THE 

QUARK-MESON-MODEL 

The quark-meson model is an SU{2)l (E) SU{2)[i invariant linear (j-modcl with chiral 
mesons (p — (cr, tt) coupled to constituent quarks q. It is an effective model for dynamical 
spontaneous chiral symmetry breaking at intermediate scales of A; < Auv, where the ultra- 
violet scale Auv ~ 1-5 GeV is determined by the validity of a hadronic representation of 
QCD. At the UV scale A^/y, the quark-meson-model is defined by the effective action 

^Auv M = J {Q^9q + gq{a + ir ■ n-f5)q + + + rricqq} (2) 

in a four-dimensional Euclidean volume with compact Euclidean time direction. The parti- 
tion function Z has the path integral representation 

Z= Jvqjvqjv(j>exp{-rj,^^), (3) 

where in the Euclidean time direction periodic and anti-periodic boundary conditions apply 
for bosons and fermions, respectively. A Gaussian approximation to the path integral fol- 
lowed by a Legendre transformation yields the one-loop effective action for the scalar fields 

0, 

m = rA,,[(/>] - IV log (rg)[0]) + ^Tr log (rg)[0]) (4) 

where rg^[0] and rP[(f)] are the inverse two-point functions for the bosonic and fermionic 
fields, evaluated at the vacuum expectation value of the mesonic field 0. Here, the boundary 
conditions of the functional integral appear in the momentum traces and we neglect con- 
tributions from mixed quark-meson-loops. In order to regularize the functional traces, we 
use the Schwinger proper time representation of the logarithms. We consider the effective 
action F in a local potential approximation (LPA), which represents the lowest order in 
the derivative expansion and incorporates fermionic as well as bosonic contributions to the 
potential density U. In this approximation, the effective field is considered to be constant 
over the entire volume. 
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The scale dependence is introduced through the infrared cutoff function /a(Tfc^), which 
regularizes the Schwinger proper time integral. A cutoff function of the form 

9 „ . ,o. 2 



Via 



22 



-{rk 



(5) 



23 



24 



. We obtain a renormalization 



satisfies the required regularization conditions 
group flow equation for the effective potential by performing a renormalization group im- 
provement and replacing the bare two-point functions by the renormalized, scale-dependent 
two-point functions, 



1 /"^ dr 
-Tr / — 

2 Jo r 

-Tr / — 



exp[-rrg),[0]] 



exp[-rrg^^ 



(6) 



In LPA, the effective action reduces to the effective potential through the relation 

Tk[<f)] = I d^x Uk{(p). (7) 



The renormalization group improved evolution equation for the effective potential in infinite 
volume is given by 



oo,r) 



dr f d^p 



m.NfeM-rip' + M2((t, vf^))] 



(2vr)4 

_exp[-r(p2 + Mlict?))] - 3exp[-r(/ + Ml{ct>'))]]k^fa{rk 



(8) 



It is necessary to choose the parameter a in such a way that the resulting integrals over the 
proper time parameter r remain finite. In infinite volume, the lowest possible integer value 
is a = 2. The diagonalization of the meson mass matrix gives the running meson masses 
which depend on the effective potential. Without explicit symmetry breaking, they are of 
the form 



Ml 
Ml 



dUk . d^Uk 



(902 

dUk 



+ 



id(l> 



2\2 



(9) 
(10) 



To derive renormalization group flow equations in the finite volume L'^ ^ at finite tempera- 
ture T, we replace the integrals over the momenta by a sum 



dpi 



(11) 



n,=— oo 
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and apply periodic boundary conditions for bosons and anti-periodic boundary conditions 
for fermions in time and space-like directions. The sums run from — oo to +00, where the 
vector n denotes (ni, ?t,2, n^^i). The Matsubara frequencies take the value ujn = 2TtnT 
for bosons and Vn = (2n + l)7rT for fermions, respectively. In the following we use the 
short-hand notation 

p^ = Erf = ^E(".4) (12) 

pI = Er-^E-.'; (13) 

i=l 1=1 

for the momenta of the fermions and bosons, respectively. In finite volume, we allow explicit 
symmetry breaking in the effective potential, which then becomes a function of two variables 
a and vf^ separately. The corresponding expression to eq. (jH)) is 

k^U,{a,n\L,T)=\^ {4iV.iV, exp [-r(.f + p| + M>, ^1] 

I n 



^exp[-r(c.,2+p| + M?^(a,vf2))] }k^fa{re). (M) 

i=l 



The Mf, z = 1, . . . , N'j, Nf = 2, are the eigenvalues of the second derivative matrix 

(^''-*^')'-||- 

of the meson potential [/^(c", tt^) with respect to the fields (f) = (a, vf). They depend only 
on the magnitude of the pion fields vf^ and are independent of the direction. We wish to 
stress the importance of this point, since otherwise the meson contributions from the flow 
equations are not compatible with the ansatz for the potential which we will introduce below. 
The second derivative matrix is given by 

/ U^^ f/^v27r(i) f/^v2 7r(2) f/^2,2 7r(3) \ 

?7,^227r« 2f/^2 + f/^2^24(7r«)2 f/^2^24 ttWtt^^) f/^2^24 tt^tt^^) 

f/,^22 7r(2) f/.2^24 7r(2)7rW 2 1/^2 + [7^2^24 (7r(2))2 U^2^2ATt(^h^^) 

Vf/.^22 7r(3) f/^2^24 7r(3)7r« f/^2^247r(3)7r(2) 2f/^2 + f/^2^24 (7r(3))2 / 

where we have suppressed the scale index k of the potential and use the abbreviations 



(16) 
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and the corresponding expressions for the higher derivatives. The eigenvalues of this matrix 
are given by 



Ml 
Ml 



1 

2 

1 
2 



2 + 4 TT^ + + \l (2 f/^2 + 4 7f2f/^2^2 

2 



16 7f2 U^^ 



2 + 4 7f^ f/^2^2 



(2 t/,^2 + 4 7im^2^2 - U^^y + 16 7r2 f/2 



;i8) 



For vanishing cross terms U^jf2, the last eigenvalue reduces to 2t/^2 + 4 7r2[/^2^2, which 
corresponds to a derivative in "radial" direction in the pion-subspace. Especially for vf^ = 0, 
the three pion modes have equal masses. We also note that the pion fields appear only in 
the combination vf^ in the eigenvalues, despite the fact that the derivative matrix contains 
terms linear in 7r^"\ The reason is that the rotational symmetry of the pion space remains 
unbroken even in the presence of explicit symmetry breaking terms in the sigma direction. 

As discussed in ||32j, we are able to perform the sum over the thermal Matsubara fre- 
quencies analytically and obtain an evolution equation which contains the Fermi-Dirac- 
distribution npiE) and the Bose-Einstein-distribution nsiE), 



k^U,{a,n ,L,T) = 2r(a + 1) L^-^ Jdk^] 

4iV,iV. 



4 ^ 

n i=l * 



with 



npiE) 



E_ 1 

er + 1 



nsiE) 



2nF{E,)) 



- 1 ' 
eT — 1 



(19) 



(20) 



where in the absence of a chemical potential for the fermions the Fermi-Dirac-distributions 
for quarks and antiquarks coincide. The effective energies are defined by 



El = k'+pl + M^{a, vf2); Ef = + pi + Mf{a, vr^). 



(21) 



The parameter a in the cutoff function is given by a = 2 in the finite volume, in order to 
ensure that we can compare with the results in infinite volume. In the limit of vanishing 
temperature T = 0, we find in d = 4 the finite volume evolution equation of the meson 
potential: 



d 

k—Ukia, TT^L) 



16 L3 



E 



f 



1 



iE,in,L))^ ^(E,(r?,L)) 



(22) 
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The polynomial ansatz for the meson potential is determined by the following idea: Since 
the current quark mass is the only source of symmetry breaking, the quark term in the 
flow equation determines the symmetry breaking terms of the potential. The constituent 
quark mass can be expanded around a finite expectation value of the mesonic fields, which 
is chosen in the direction of the field a, 

= /[(a + m^f + TT^] = / [(a + ao - (Tq + m^f + vr^] 

= /[((To + + 2m,((T - ao) + (a' + vf^ - a^)]. (23) 

We have rescaled rric by a factor g for convenience, so that the physical current quark mass 
is given by grric. From this expression, we read off that the contributions to the potential 
from the fermionic terms in the flow equations can all be expressed in terms of powers of the 
combinations (cr^ + vf^ — a^) for the symmetric part and {a — ctq) for the symmetry breaking 
parts. Therefore, we make for the meson potential the ansatz 

U,{a,n') = Y, E ci,,{kKa-aoy{a' + n'-a',y (24) 

1=0 j=0 

The flow equations for the coefficients in this potential are derived in the appendix. 

Incorporating the explicit breaking of the chiral symmetry into the potential and the 
flow from the start has several advantages. The polynomial expansion above evolves au- 
tomatically from a potential with small symmetry breaking peaked around (a) ^ to a 
potential with large symmetry breaking peaked at a value (a) ~ Without explicit sym- 
metry breaking, the polynomial expansion in 0^ has to be changed from a parametrization 



in terms of powers of 0^ to (0^ — 0g) at the chiral symmetry breaking scale 22 1. 

As is well known, a linear symmetry breaking term remains unchanged in the renor- 
malization group flow js^. Therefore the usual strategy is to evolve the potential without 
a symmetry breaking term. Explicit symmetry breaking is then taken into account after 



29, In an infinite 



the quantum fluctuations have been integrated out on all scales |2^ 
volume and for small quark masses, this is perfectly acceptable and will yield the correct 
results. In a finite volume, however, the situation is different. Since chiral symmetry is not 
spontaneously broken, explicit symmetry breaking has to be included on all scales in the 
renormalization group flow to obtain a nonzero value for the order parameter. Otherwise, 
divergences from massless Goldstone bosons would restore the symmetry. In this context. 
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we would like to point out that even in the absence of a symmetry breaking term, the pion 
decay constant does not remain zero on all renormalization scales k. On some intermediate 
scale below the chiral symmetry breaking scale, k < k^, where the quantum fluctuations are 
only partially integrated out, it acquires a nonzero expectation value, and chiral symmetry is 
spontaneously broken. However, the emergence of exactly massless Goldstone bosons dom- 
inates the infrared evolution of the potential and counteracts the formation of a symmetry 
breaking condensate. 

When the potential is expanded in a polynomial in a theory with exactly massless Gold- 
stone bosons, divergences appear in the flow equations for the coefficients of operators of 
mass dimension higher than four [2^. As an added benefit of including explicit symmetry 
breaking, the presence of a finite pion mass regulates these IR divergences. 



IV. NUMERICAL EVALUATION 



We have solved the RG fiow equations numerically and present the results for the volume 
dependence of the pion mass and the pion decay constant in the following section. For the 
numerical evaluation, we have used the polynomial ansatz for the effective potential given 
in eq. ()24|) . and expanded up to fourth order in the fields: 

Uk{<7, n^) = aoo{k) + aoi(fc)(a^ + vf^ - o-q) + (^02{k){cr'^ + tt^ - a^ f 

+aio{k){a - (To) + a2o{k){a - aof + a3o{k){cr - aof + a4o{k){cr - (Jq)^ 
+au{k){a - ao){a^ + 71^ - al) (25) 

Here, we first discuss our choice of model parameters at the UV scale, and some details of 
the numerical evaluation. 

The UV scale itself is determined from physical considerations as the scale below which 
a description of QCD with hadronic degrees of freedom is appropriate. Here, we choose 
Ajjv = 1-5 GeV. At the ultraviolet scale A^y, the free parameters of the quark- meson- 
model are the meson mass muv, the four-meson-coupling Xuv, and the current quark mass 
grric, which controls the degree of explicit symmetry breaking. The Yukawa coupling g does 
not evolve in the present approximation j^, choose g = 3.26, which leads to a 

reasonable constituent quark mass of Mg = glf^^ + rric) ~ 310 MeV for physical values for 
the pion decay constant f^, = 93 MeV and the current quark mass grric = 7 MeV. 
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Kuv [MeV] muv [MeV] Xuv 5"^c [MeV] [MeV] [MeV] 

1500 779.0 60 2.10 90.38 100.8 

1500 747.7 60 9.85 96.91 200.1 

1500 698.0 60 25.70 105.30 300.2 



TABLE I: Values for the parameters at the UV-scale used in the numerical evaluation. The 
parameters are determined in infinite volume by fitting to a particular pion mass and the corre- 
sponding value of the pion decay constant, taken from chiral perturbation theory. Note that in our 
notation, the physical current quark mass corresponds to grric- 

In table n we summarize the three parameter sets which we used in obtaining our results 
for pion masses of 100, 200 and 300 MeV. We determine these UV parameters by fitting to a 
particular value for the pion mass m^(oo) and to the corresponding value for the pion decay 
constant firioo) in infinite volume. We then evolve the RG equations with these parameters 
to predict the volume dependence of f-jriL) and mT^{L). 

For any value of the pion mass, the corresponding value of the pion decay constant is 

n 

taken from chiral perturbation theory The pion mass is mainly controlled by the value 
of the current quark mass, which parametrizes the symmetry breaking. The current quark 
mass varies from approximately 2 MeV for a pion mass of 100 MeV to about 10 MeV for 
m^r = 200 MeV, it has to be increased to approximately 25 MeV for m.,^ = 300 MeV. To 
achieve the correct corresponding values for the pion decay constant, the meson mass at the 
UV scale has to be decreased from approximately muv = 780 MeV to muv = 700 MeV, 
while the pion mass increases from 100 to 300 MeV. The four- meson-coupling Xuv is fixed. 
We have checked that our results are to a very large degree independent of the particular 
choice of UV parameters: Different sets of parameters leading to the same values of the 
low-energy constants in the infinite volume give the same volume dependence. 

Although it facilitates the comparison to chiral perturbation theory, it is not necessary as 
a matter of principle to use the chiral perturbation theory result for the mass dependence of 
the pion decay constant. However, as has been found for infinite volume, in order to correctly 
describe the behavior of the pion decay constant as a function of a single symmetry breaking 
parameter, it is necessary to go beyond the approximation of a constant expectation value for 
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the meson field, which we used in this paper, and to include wave function renormalizations in 
the RG flow 23]. This makes it possible to recover the correct prefactors of chiral logarithms 
in the framework of the renormalization group. Such an approach is more powerful than the 
present one, since in addition to the volume dependence, it predicts the dependence of 
and fjr on the symmetry breaking parameter rric- We stress that even in such an approach, 
it remains necessary to fit the parameters at the UV scale to the correct values of the low 
energy constants. Thus, for example the value of the pion decay constant in the chiral limit 
is not a prediction of the model, but a necessary input to constrain its parameters. The 
full set of RG-equations including the wave function renormalization and coupling constant 
renormalization equations would reduce the input parameters to the four-fermion coupling 
and the current quark mass at the UV-scale. In connection with the symmetry breaking 
ansatz eq. (j211), these equations are more complicated and have not been worked out yet. 

A limit on the possible values of the current quark mass is given by the requirement that 
all masses, in particular the sigma-mass, remain substantially smaller than the ultraviolet 
cutoff Auv ~ 1500 MeV of the model. For a pion mass of 771^^ = 300 MeV, we find 
^ 800 MeV. 

With regard to the UV-cutoff, we find only a slight dependence of our results for reason- 
ably large volumes. When we change the cutoff from A^y = 1500 MeV to Auv = 1100 MeV, 
our results for the relative shift i?[m^(L)] of the pion mass in the finite volume change little. 
The change in the pion mass from a variation of the cutoff is of the order of less than 1% for 
L > 2 fm, and approximately 6% at L = 1 fm for the largest pion mass we considered here, 
rTT-TT = 300 MeV. For smaller pion mass, the dependence on the UV cutoff becomes weaker, 
for = 100 MeV it is negligible on the scale of our results. This can be understood, since 
a higher degree of explicit symmetry breaking leads to more massive particles for which a 
smaller value for the UV momentum cutoff becomes more relevant. 

The sums over the momentum modes in the flow equations cannot be performed ana- 
lytically. For a numerical evaluation of the flow equations, these sums must be truncated 
at a maximal mode number Amax = max|fi| which deflnes the cutoff momentum mode 
Pmax = x^™ax- This numerical truncation should not introduce an additional UV cutoff in 
the model, and therefore we require that 

^AT^ax > Auv- (26) 
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FIG. 1: The figures show the dependence of the results for the pion mass in finite volume on the 
truncation A'^max of the momentum sum in the numerical evaluation of the flow equations. Plotted 
is R[mTj{L)], the relative deviation of the finite volume pion mass from the value in infinite volume, 
as a function of the box size L, for different values of iVmax- Results for two different values of the 
pion mass are shown in the two panels. Ideally, for large volumes, the relative deviation from the 
infinite volume pion mass should approach zero. The truncation has a relatively larger effect for 
heavier pion masses. 

Since we use a "soft" cutoff function it is necessary to really satisfy the above equation with 
a safe margin. For the volumes with L > 1 fm considered here, we have used A^max = 40. 

In figure ^ the dependence of the relative difference of the pion mass in finite volume 
from its value in the infinite volume limit, i?[m^(L)], is shown as a function of L for different 
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values of the maximal mode number A^max- The relative mass difference depends mainly on 
m-j^L and drops exponentially for large values of this dimensionless variable. Thus, for any 
given value of L, the value of R[m.,t{Ly\ will be smaller for a heavier pion. Comparing the 
two panels in figure ^ we see that although the absolute values are smaller, the relative 
error due to the finite number of momentum modes in the evaluation is larger for a heavier 
pion. The reason is the increasing importance of the non-zero momentum modes when 
the pion mass becomes larger at fixed box size. Following the argument from ^| outlined 
above, if ^ L, then the partition function is dominated by the zero modes and effects 
from finite momentum modes present small corrections. When the pion mass is increased, 
the importance of the finite momentum modes grows and the number of modes has to be 
increased to obtain results with the same level of accuracy. This argument can be presented 
in a more formal way. The momentum sums contributing to the flow equations are of the 
form 

rii,n2,n3 ' tt ' i^-^ 

where "n? = n\ + n\ + n\. For small values of the renormalization scale fc, the sum is 
dominated by the zero mode term 



If the box is sufficiently small, all terms with non-zero momentum are suppressed by l/i^^, 
which acts as a large regulator. As we increase the size of the box, so that 1/L ~ itLt^^ the 
contributions of the nonzero momentum modes are of the same size as the zero mode term. 
Therefore, effects from a truncation of the momentum sums should be expected to appear 
already at a smaller box size for large pion masses, which is exactly what we observe. 



V. RESULTS 



The results of the RG-flow equations for the evolution with the infrared cutoff scale 
k give a picture of chiral symmetry breaking which reflects the formation of the quark 
condensate for higher momenta and the effects of pion fluctuations at low scales. Figure El 
shows the masses of the pion and sigma, and the pion decay constant as a function of 
the renormalization scale k, for a value of = 100 MeV, in the infinite volume limit. 



15 




k [MeV] k [MeV] 

FIG. 2: Masses of the mesonic degrees of freedom and the pion decay constant as a function of 
the renormaUzation scale k in infinite volume. The chiral symmetry breaking scale can be clearly 
identified as the scale at which the mass of the heaviest meson (the a) has a minimum. For this 
figure, we have chosen m.jr{oo) = 100 MeV and /7r(oo) = 90.4 MeV. 




1/L [fm"''] 1/L [fm"'] 



FIG. 3: Masses of the mesonic degrees of freedom and the pion decay constant as a function 
of the inverse box size l/L. The results are obtained by completely integrating out all quantum 
fluctuations {k — > 0) at fixed L. As soon as k < 1/L, the box size becomes the controlling scale, 
and in the limit A; it is the only scale that remains. As for the preceding figures, we show the 
results with = 100 MeV and /^r = 90.4 MeV for A; ^ and L ^ oo. 

Starting at the UV scale Auv and proceeding towards smaller values of A;, we observe that 
the pion decay constant /^r grows rapidly around the chiral symmetry breaking scale k^^sB ~ 
800 MeV, begins to flatten between 600 — 400 MeV, and becomes almost completely flat 
below 300 MeV. Generally, massive degrees of freedom decouple from the renormalization 
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group flow at a momentum scale given by the value of their mass m, i.e. they do not 
contribute to the renormalization for k < m. This can be seen clearly in the flow of f^^. As 
soon as the renormalization scale is of the order of the constituent quark mass (approximately 
300 MeV), the quarks are no longer dynamical degrees of freedom and /tt becomes essentially 
constant. The RG flow of the mass of the heaviest meson, the sigma, is in several respects 
very similar to the flow of Jt^. Its slope is also initially large at the chiral symmetry breaking 
scale and starts to decrease between 600 — 400 MeV as well. The value of the sigma- 
mass reaches a maximum at k slightly above 300 MeV. Its decrease below this scale is 
due to the light pion with a mass of 100 MeV, which remains in the evolution as the only 
dynamical degree of freedom. When the pion mass is increased, the drop in the sigma- 
mass below the scale set by the constituent quark mass becomes much less pronounced. 
For ^ rriq ^ 300 MeV, m^{k) is essentially a flat function of k after it has reached its 
maximum. 

In flnite volume, a similar behavior is visible. In flgurelHl the meson masses and the pion 
decay constant are shown as a function of the scale 1/L set by the flnite volume. In these 
results, all quantum fluctuations are integrated out completely, which removes the scale k. 
Now let us consider a flnite value of k, where the quantum fluctuations are only partially 
integrated out. The scale 1/L introduced by the flnite volume is in competition with the 
renormalization scale k. As soon as k drops below vr/L, the renormalization scale no longer 
controls the renormalization flow. We can interpret the results shown in flg 01 roughly as 
an instant picture of the k-flow arrested at a scale k = n/L. However, this correspondence 
is not one-to-one: while the cutoff k affects both bosonic and fermionic flelds in the same 
way, this is not true for 1/L. Since there are no zero modes for flelds with anti-periodic 
boundary conditions, the fermionic flelds are more strongly affected by this cutoff than the 
mesons. For the mesons, the scale ^ imposes only a minimum value for the smallest non-zero 
momentum mode. For the fermionic flelds, on the other hand, the lowest momentum mode 
can effectively "freeze" the quark flelds already above the constituent quark mass 
scale and no condensation of quarks takes place. For very small volumes 1/L > 0.5 fm^^, 
the suppression of quark condensation by the large cutoff becomes the dominating effect. 
The chiral symmetry is approximately restored. A more subtle effect can also be seen in the 
behavior of the sigma-mass. While the sigma-mass has a maximum in the fc-flow at a value 
of ^ 600 MeV, from which it drops to m„ ^ 500 MeV due to the pion fluctuations. 
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there is no corresponding maximum in the 1/L-dependence. With decreasing 1/L the pions 
hardly feel the constraints of the finite volume, but the momentum scale at which the quarks 
decouple consistently increases. Therefore the pion contributions at low momenta have a 
greater effect in the RG flow, since the /c-region increases in which they are the only relevant 
degrees of freedom. Through this effect, the pions also contribute toward the restoration of 
chiral symmetry for small volumes. 

The volume dependence of the low energy observables is the main results of this paper. 
In table m we give the values for i?[m^(L)], cf. eq. (0), the relative difference of the pion 
mass m.^{L) in finite volume from m.^{oo), its value in infinite volume, for three volume sizes 
L = 2.0, 2.5, 3.0 fm and three pion masses m^(oo) = 100, 200, 300 MeV. 

In figures m and we show the relative change of the pion mass mT^{L) and the pion 
decay constant f.j^{L) as a function of the size L of the three-dimensional volume. We plot 
the results for the relative differences on a logarithmic scale for three different values of the 
pion mass m^(cx)) = 100, 200, 300 MeV. 

Let us first discuss the plots for the pion mass in fig. |31 The relative change of the pion 
mass decreases with the volume size L and the pion mass m,r(oo). Fig. 0] also contains the 
results of chiral perturbation theory from the exact one-loop calculation in finite volume 
[lol | and also the "best estimate" from f9|. As discussed in section m the difference between 
the results of Liischer's formula for the mass shift which uses as input Trvr-scattering in the 
one-do) and three- (nn/o) loop order is used as a correction to the exact one-loop mass shift 
from |lO| . 

Our RG results have the same slope as those from chiral perturbation theory, but are 
consistently above chiral perturbation theory, even with corrections to three loops. In gen- 
eral, the RG calculation gives values for Rlm-^lL)] which are about a factor 1.5 to 2.0 larger 
than the values from chPT, as can be seen in table ITTl 

The difference between the RG result and the loop expansion decreases with higher order 
in loops. The RG result is closer to the calculation with nn/o-Liischer formula than to the 
one-loop chPT calculation. 

For large volumes, the pion mass should drop as exp{—mT^L). Therefore, we expect that 
the slope of the RG result is the same as that from chiral perturbation theory, which is 
indeed the case. The calculation with the flow equation can be continued to smaller box 
sizes, as long as the momenta constrained by the box are below the cutoff Auv The slope 
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R[m^{L)] 



L [fm] 


m^(oo) [MeV] 




RG 




IL chPT 




+(nnlo-lo) 


AR 




2.0 


100 


1.01293 


26.6 X 10" 


-2 


8.74 X 10" 


2 


11.6 X 10-2 


15.0 X 10" 


-2 




200 


2.02586 


5.38 X 10" 


-2 


2.00 X 10" 


2 


3.31 X 10^2 


2.07 X 10" 


-2 




300 


3.03879 


1.70 X 10- 


-2 


0.56 X 10- 


2 


1.12 X 10-2 


0.58 X 10- 


-2 


2.5 


100 


1.26616 


10.37 X 10" 


-2 


3.85 X 10" 


2 


4.97 X 10-2 


5.40 X 10" 


-2 




200 


2.53233 


1.95 X 10" 


-2 


0.73 X 10- 


2 


1.17 X 10-2 


0.78 X 10- 


-2 




300 


3.79849 


5.31 X 10" 


-3 


1.65 X 10^ 


3 


3.27 X 10-3 


2.04 X 10" 


-3 


3.0 


100 


1.5194 


4.94 X 10" 


-2 


1.91 X 10" 


2 


2.41 X 10-2 


2.53 X 10" 


-3 




200 


3.03879 


7.85 X 10' 


-3 


2.95 X 10" 


3 


4.65 X 10-3 


3.20 X 10" 


-3 




300 


4.55819 


1.76 X 10" 


-3 


0.54 X 10" 


3 


1.05 X 10-3 


0.71 X 10- 


-3 



TABLE II: Values for R[m-,^{L)], cf. eq. the relative shift of the pion mass in finite volumes 



of L = 2.0, 2.5, 3.0 fm, compared to the value in infinite volume, for pion masses of 



oo 



100, 200, 300 MeV. We compare our RG calculation to the exact one-loop chPT results of [llj for 
a finite volume (IL chPT), and the exact one-loop calculation with corrections in three- loop order 

n 

obtained with chPT using Liischer's formula 9] (IL chPT -|- {nnlo-lo)). In the last column, the 
difference Ai2 between the RG result and the three-loop corrected chPT result is given. 

of the RG result at small box size (L < 1 fm) is approximately given by the meson mass 
TTirj = m^, cf. Fig. 121 after the transition between the regime dominated by chiral symmetry 
breaking and the one with restored chiral symmetry has taken place. In this other region 
the chiral expansion can no longer be considered reliable and is therefore no longer applied. 

While the relative difference between the exact one-loop result and the RG results remains 
approximately constant for different pion masses, the relative difference between the RG 
results and the results of Colangelo and Diirr jj^l decreases when the mass of the pion is 
increased. We have checked that the difference between our RG results and the chiral 
perturbation theory results from Colangelo and Diirr are consistent with the error estimate 
of Liischer's approximation formula. We find for all pion masses that the differences decrease 
exponentially according to exp(— C m^L), with C a positive constant of order 1. 

Next we discuss the results for the volume dependence of the pion decay constant shown 
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FIG. 4: Volume dependence of the pion mass. We plot the relative shift of the pion mass from its 
infinite volume limit = {mT^{L) — m7r(oo))/m7r(cxD) as a function of the^ize of the volume 

L. For comparison, we also plot the results from chPT calculations taken from 9]. The values for 
the pion mass in infinite volume are given in the figure, note the different scales on the axes for 
different m7r(oo). 
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FIG. 5: Volume dependence of the pion decay constant. We plot the relative shift from the infinite 
volume limit Rf^{L) = (/^(oo) — fT^{L)) / fT^{co) as a function of the box size L for different values 
of the pion mass. For comparison, we show the chPT results taken from Q]. The values for 771,7^(00) 
are given in the panels. Note the different scales on the axes for different 711,7^(00). 
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in fig. 13 As for tlie pion mass, we compare witli tlie cliiral perturbation tlieory results 
jl^ and j^. Note that in this case, we define the relative difference with opposite sign 
as R[fT,{L)] = (/7r(oo) — /^(L))//7r(oo), because the pion decay constant is smaller for a 
finite volume, in contrast to the pion mass. We refer to fig. Elfor an illustration of the global 
behavior of fn{L). When the size of the volume is decreased, the pion decay constant at first 
drops only slowly, and then sharply at the scale associated with chiral symmetry breaking. 
For very small volumes, the pion decay constant vanishes and chiral symmetry is effectively 
restored. This is true for L — » regardless of the value of the pion mass. Therefore, the 
largest possible relative shift for L ^ is /^[/^^(O)] = 1, when the order parameter vanishes 
completely. 

The plots for the volume dependence of /,r illustrate in which region chiral perturbation 
theory remains valid. Since spontaneously broken chiral symmetry and a sufficiently large 

are required, chPT does not describe the transition to the region with effectively restored 
chiral symmetry. Therefore, chiral perturbation theory results are available only for volumes 
with L > 2 fm. Already for L = 2 fm and a pion mass = 100 MeV, the /7r-shift is almost 
40%. The RG method remains valid in both regions, down to very small volumes. 

We compare with the chiral perturbation theory results obtained by using an approxima- 
tion similar to Liischer's formula for the pion mass shift, which was derived in For the 
pion decay constant, the input needed to calculate the finite volume shift is an amplitude 
involving the axial current in infinite volume. So far only two loops in chiral perturbation 
theory {nlo) are known. The chPT n/o-corrections to the exact one- loop result increase the 
shift towards larger values. As in the case of the pion mass, the RG results give a slightly 
larger finite volume shift than the chPT results. Again the results from RG and chPT con- 
verge for larger values of the pion mass (note the different scales on both axes in the plots 
for different values of 771^^1). 

For the volume dependence of the sigma-mass, we refer back to fig. El where we show the 
overall dependence of rrifj on 1/L. It is interesting to relate the variation of mcr{l/L) to the 
0++ phase shift of Trvr-scattering, which will be done in a separate work. 
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VI. CONCLUSIONS 

We have presented a new approach to the quark-meson-model employing the renormal- 
ization group method in a finite volume within the framework of the Schwinger proper time 
formalism. Central to any such approach is the inclusion of explicit chiral symmetry break- 
ing. Since chiral symmetry is not broken spontaneously in a finite volume, it is necessary to 
introduce a finite current quark mass. In this paper, we have evolved the effective potential 
with additional symmetry-breaking terms. The form of these terms is constrained by the 
quark contributions to the renormalization group fiow, which introduce the explicit chiral 
symmetry breaking. 

By solving the resulting renormalization group fiow equations numerically, we have ob- 
tained results for the volume dependence of the meson masses, in particular the pion mass, 
and the pion decay constant, the order parameter of chiral symmetry breaking. 

Our results show consistently a larger finite volume mass shift for the pion than has been 
obtained in chiral perturbation theory including up to three loops. The differences between 
the chiral perturbation theory results which make use of the Liischer formula and the RG 
results obtained in the present paper are consistent with the error estimate for Liischer's 
approximation. As one expects, the difference is largest for small values of itIt^L. We have 
checked that this difference decreases exponentially with an increase in this dimensionless 
quantity. As shown in figs. El our results and those obtained in chiral perturbation theory 
with Liischer's formula [a, Isl converge for large current quark masses. We note that the 
ratio of the results from chPT and RG does not depend on L, even down to L = 1.5 fm. 

Compared to the present numerical approach, chiral perturbation theory has the advan- 
tage that it is possible to obtain analytical expressions for the finite volume mass shift. This 
makes comparisons to lattice results simpler. On the other hand, current lattice volumes 
and lattice pion masses are at the edge of a reliable chiral perturbation theory calculation. 
In contrast, the RG method remains valid for large current quark masses as well as small 
volumes. 

The main uncertainty of the RG method comes from its dependence on the UV cutoff 
scale A[/y for large meson masses. The system becomes sensitive to Kuv for large explicit 
symmetry breaking, because the mass of the sigma as the heaviest particle approaches the 
UV cutoff. For a pion mass of ttIt^ = 300 MeV, the sigma mass is ma ~ 800 MeV. In this 
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case, a cutoff variation between Kjjy = 1500 MeV and 1100 MeV, changes the pion mass for 
a volume with L = 1 fm by approximately 6%, and by less than 1% for L > 2 fm. Within 
this uncertainty, our results agree with those of chPT for = 300 MeV. In contrast, for 
m^r = 100 MeV the cutoff dependence is so weak that it is not noticeable on the scale of 
the results. The RG and chPT results do not agree within this uncertainty, cf. Table ITTl for 
= 100 MeV and = 200 MeV. 

The dependence of our results on the choice of model parameters at the UV scale is much 
weaker than that on the cutoff. By fitting to the values of the low-energy observables m„ 
and /tt in infinite volume, we achieve a very high degree of independence on the particular 
choice of UV parameters. 

We expect that the inclusion of wave function renormalizations in the finite volume 
renormalization group flow should make it possible to describe the low energy constants as 
a function of a single symmetry breaking parameter, the quark mass m^. It has already 
been observed that this is the case in infinite volume, where the behavior of m.^ and f^r 



as a function of rric is described correctly when all other UV parameters remain fixed [27]. 
Systematic errors introduced by the simplification used in this paper to adjust both m7r(oo) 
and fnioo) in infinite volume could then be estimated by a comparison of the calculations. 

Although we have not yet investigated these questions in depth, it appears that the 
present approach, which treats the pion fields and the sigma explicitly and as individual 
degrees of freedom, improves the convergence j^] of the polynomial expansion of the effective 
potential. This is due to the finite mass acquired by the Goldstone bosons, which becomes 
relevant for A; — > 0. 

The RG approach shows in a transparent way the relevance of the momentum scale 
introduced by the finite volume for the quantum fluctuations. 




Acknowledgments 



H.J. P. would like to thank Prof. Leutwyler and Dr. Diirr for instructive discussions. J.B. 
would like to thank the GSI for financial support. 



24 



APPENDIX A: FLOW EQUATIONS 

We derive flow equations for the coeflicients of the potential by inserting eq. ()24|) into the 
flow equation for the potential and comparing coefficients of both sides. We define expansion 
coefficients for the flow equations: 



"dk'^y z!j! [da J [dn^J "dk'^' 



(Al) 



Then the flow equations for the coefficients Oj,- in the ansatz for the potential are given by 



5 \ f d \ f da 



+ 1) (-2^oA:^) (1 - S.iN.-.),) (A2) 

with the additional condition that (1 + 2j) < N^j. 

In order for o" = ctq to actually correspond to a minimum of the potential, it must satisfy 
^ f^fclo-=(To ~ ^- -^^^ coefficients aio and agi, this translates into the condition 

aio + 2aoio-o = 0. (A3) 

Due to this condition, only two of the variables in the set {aio, o-oi, ctq} are independent, and 
the third one can be expressed in terms of the other two. Likewise, if we take the derivative 
of the equation ()A3|) with respect to the renormalization scale k, we get an equation which 
relates the flow of these three variables: 

d f d \ d 

^'dk"'^^ I dk^^ J ^^^^'dk^^^ ^ ^' ^^^^ 



We can use this equation to replace the flow equation for aio in the above set eq. ()A2|) . 
It is desirable to eliminate aio, since ctq and aoi both correspond to the observables we 
wish to obtain, namely the pion decay constant and the pion mass. In addition, with this 
replacement the system of differential equations can also be solved more easily. 

From the general expression for the flow equations eq. ()A2jl . we flnd the particular equa- 
tions governing aio and aoi: 

d d f d \ f d \ 

^df' ) = - y'dk''' ) - ^2-°^^"° ) ^^^^ 

d 8 f . d \ f . d 



f^Q^Uk 1 = ^^«oi - an ( ^^^0 1 - 2ao2 ( 2aok—ao ] . (A6) 
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The flow equation for aiQ contains on tlie LHS only terms that are proportional to the 
symmetry-breaking current quark mass rric. Because of this, aio does not evolve in the 
chiral limit rric — > 0. If it is initially zero at the UV scale, it remains zero on all scales. In 
this case, the condition ()A4jl forces the coefficient aoi to vanish as soon as (Tq acquires a 
finite expectation value. This corresponds to the appearance of exactly massless Goldstone 
bosons in case of spontaneous symmetry breaking, in accordance with our expectations for 
the chiral limit. 

In order to derive a fiow equation for the minimum of the potential ctq, we can combine 
the two equations and use eq. ()A4j) to eliminate the /c-derivatives of aio and aoi: 



d..V° i.a 



d \ 

A;^o-o 1 (2a2o + 2aoi + 4aiiao + 8ao2cro) . (A7) 



From the expressions for the meson masses, evaluated at the minimum of the potential, 
it can be seen that the expression in brackets, which multiplies the /c-derivative of (Tq, is 
up to a constant factor the square of the a-mass, M^. Therefore, this equation is always 
well-conditioned. The only exception is at the chiral symmetry breaking scale, where 
drops sharply, if the explicit symmetry breaking is very small. For reasonably large pion 
masses, this is not a problem. 

APPENDIX B: ANSATZ FOR THE POTENTIAL 

We use the fiow equations in infinite volume to motivate the ansatz for the potential with 
explicit symmetry breaking. Neglecting the mesonic contributions to the fiow equations, we 
are left with the terms arising from the fermions: 

k-u-(.y) = -^p^. (Bi) 

The constituent quark mass Mq contains through the current quark mass nic the only ex- 
plicitly symmetry breaking term in the fiow equation. As shown in eq. ()23|). by expanding 
around the minimum of the potential, the quark mass can be written as 

Ml = g\ao + m,f + 2m,{a - (Tq) + (a^ + yr^ - a^)], (B2) 

where nic is rescaled by a factor of g for convenience. When we expand the denominator 
in the fiow equation in the deviation of the fields from the vacuum expectation value, the 
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result contains only those terms we postulated in our ansatz for the potential: 



327r2 



ANfNc 



X 



I I -;2 2\2 

+ ((7 +7r -Gq) 



9' 



2 



+ ((7 - (7o) 
+ ((7 - (7o)' 



—2mr 



Ami 



2 



A;2 + (yf2(o'o + mc)2 

2 \ 2 



5 



^ ^ A;2 + 52(^Q+^^)2 



+((7-(7o)((7" + r-ao') 2 



—2mr 



+ 



(B3) 



All remaining terms are of higher order in {a — o"o) or (cr^ + tt^ — cxo) or any combination 
thereof. Only the terms in (cr^ + tt^ — Uq) = (0^ — ctq) remain in the chiral limit rric 0. 
These terms respect the chiral symmetry and the potential reduces to the symmetric form. 
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